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ABSTRACT 

Star formation efficiency controlled by the protostellar outflow in a single 
cloud core is investigated by three-dimensional resistive MHD simulations. Start- 
ing from the prestellar cloud core, the star formation process is calculated until 
the end of the main accretion phase. In the calculations, the mass of the prestel- 
lar cloud is parameterized. During the star formation, the protostellar outflow 
is driven by the circumstellar disk. The outflow extends also in the transverse 
direction until its width becomes comparable to the initial cloud scale, and thus, 
the outflow has a wide opening angle of ^ 40°. As a result, the protostellar 
outflow sweeps up a large fraction of the infalling material and ejects it into the 
interstellar space. The outflow can eject at most over half of the host cloud mass, 
signiflcantly decreasing star formation efficiency. The outflow power is stronger 
in clouds with a greater initial mass. Thus, the protostellar outflow effectively 
suppresses star formation efficiency in a massive cloud. The outflow weakens sig- 
nificantly and disappears in several free-fall timescales of the initial cloud after 
the cloud begins to collapse. The natal prestellar core influences the lifetime and 
size of the outflow. At the end of the main accretion phase, a massive circum- 
stellar disk comparable in mass to the protostar remains. Calculations show that 
typically, ~ 30% of the initial cloud mass is converted into the protostar and 
~ 20% remains in the circumstellar disk, while ~ 40% is ejected into the inter- 
stellar space by the protostellar outflow. Therefore, a single cloud core typically 
has a star formation efficiency of 30 — 50%. 
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Introduction 



Many protostellar outflows are observed in star-forming regions. Tliese outflow s are be- 



lieved to be universally driven by a protostar and play a critical role in star format ion (lArce et al. 



20071 : iBally et al.l 120071 ). The protostellar outflow may determine star formation efficiency, 



especially in low-mass star formation process. The similarity between the core mass func- 



(Motte et al. 


1998: 


Andre et al. 


2009) 



Recent observa- 
tions have shown that, in a single prestellar core, star format i on effi c iency is limited to 
e = Mstar/Meore ~ 20 - 50% jAudre et allboiol : iKonvves et al.lboioh . iMatzner fc McKee 



(120001 ) argued that star formation efficiency might be limited by the protostellar outflow. 
With a simple analytical approach, they showed that a wide-opening- angle outflow sweeps 
up the gas in the infalling envelope and ejects it into the interstellar space. They concluded 
that the protostellar outflow can limit star formation efficiency to e ~ 30 — 50%. In addi- 
tion, feedback from the pro tostellar outflow may maintain interstellar turbulence and affect 
subsequent star formation (iNakamura &: Lil 120071 ). Thus, the protostellar outflow is crucial 
in (low-mass) star formation. However, it is quite difficult to model the outflow using only 
an analytical approach because the outflow driving depends on various conditions such as 
the infalling mass rate onto the circumstellar disk, the size of the circumstellar disk, the 
conflguration and the strength of the magnetic field around the protostar. Thus, numeri- 
cal simulation is necessary to understand the driving and evolution of the outflow during 
star formation and to determine star formation efficiency, as controlled by the protostellar 
outflow, in more detail. 

Some authors have used a spherically symmet ric calculation to investigate the evolu- 



tion of a prestellar cloud until protostar formation (iLarsonl Il969l : IWinkler fc NewmanI Il980 



Masunaga et al. lll998l : lMasunaga fc Inutsuka 1120001 ). Now, we can directly calculate the star 
formation fr om the prestellar cloud stage until protostar formation with a multi-dimensiona l 



calculation ( lBatelll998l . l2010l : lTomisakall2002l : lBanerjee fc Pudritz 



opening-angle outflow (ITomisaka 



Tomida et al.ll2010l : ICommercon et al. 



2002 



Machida et al. 

"1 



20101: 



2006 



2005bl: 



Biirzle et al 



Machida et al.ll2007h . 



Before protostar forma tion, the first (adiabatic) core forms (Larson 1969h and drives a wide- 



lennebelle fc Teyssierll2008b 



201l[ ). After protostar form ation. 



the first core evo 



1998 



2010 



2011 



Machida et al. 



v es directly into the c ircumstellar disk after the protostar formation ( Bate 



20111 ) 



Walch et al.l l2009al : Machida et al.l l2010al : Machida fc Matsumotol 12010 



Since the first core transitions smoothly to the circumstellar disk, the 
outfiow driven by the first core before protostar formation is driven by the circumstellar disk 
after protostar formation without a transient disappearance. In summary, a wide-opening- 
angle outfiow appears before protostar formation and continues to be driven after protosta r 
formation until the end of the main accretion phase ( iMachida et al.ll2008bl . l2009al |2010b|). 
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This type of outflow i s believed to correspond to tfie molecular outflow frequently observed in 



star- f orming; regions (iBontemps et al.l 119961 : iWu et al. I I2004J : lArce et al.l 120101 : iCurtis et al. 



20101 ). Recen tly, a bipolar molecu lar outflow was observed around the candidate first core 
Per-Bolo 58 (IDunham et al.ll201in. I n addition, some fi r st core candidates were reported 



by several authors (jChen et al. 



2010l : IChen fc Arcell2010l : Enoch et al.ll2010[ ). Thus, to es 



timate the mass ejection rate from the host cloud or star formation efficiency, we need to 
calculate the evolution of the protostellar outflow from the prestellar core stage because a 
wide-opening-outflow appears before protostar formation. 

In addition to the wide-opening - angle outflow, well coUimated hig h-velocity jets are 
observed in the star forming region (lArce et al.l 120071 : iBally et al.l 120071 ). Such a jet also 
appe ars during star formation process, and is driven in very close p roximity to the proto- 
star jTomisakalboO^ : iBaneriee fc Pudritzlbood : iMachida et al.lboOSal lbl). The driving source 
and its spatial scale of each flow (wide-opening-angle outflow and well coUimated jet) are 
considerably different: the driving source for the wide-opening-angle flow is the first core or 
the circumstellar disk, which has a size of 1 AU, whereas the driving source for the well 
coUimated jet is the p rotostar and the disk near the protostar in the region of <^ 0.1 AU 
( jMachida et al.ll2008bl ). In addition, the mass ejection rate of the wide-opening-angle flow is 
much higher than that of the well coUimated jet. Because the jet has a good coUimation, it 
cannot accumulate sufficient mass in th e infalling envelope to reduce star formation efficiency 
( iTomisakal l2002t iMachida et al.ll2008b( ). Note that the difference in coUimation between the 
wide-opening-angle outflow and well coUimated jet is caused by the driving mechanism ari d 
the configuration of the magnetic field lines around the driving source (IMachida et al.ll2008bl ). 



A statistical study of molecular outflows observed in various star forming regions also 
indicates that the (wide-opening-angle) molecular outflows have a sufficient en ergy or mo- 
mentum to control the star formation (IBontemps et al.lll996l : lArce et al.ll2010l ). Therefore, 
the wide-opening-angle outflow and not the weU coUimated jet is expected to be important 
in determining star formation efficiency. Thus, to estimate the mass ejection rate (or star 
formation efficiency), we must resolve the driving source of the wide-opening-angle outflow, 
whereas we do not necessarily need to resolve the protostar itself, which drives the well col- 
limated jet. Since previous studies spatiaUy resolved the protostar, they could not directly 
estimate the resulting raass ejection rate by the protostellar outflow and star formation ef- 
flciency ( iTomisakal l2002l : iBanerjee fc Pudritzl l2006l : IMachida et al.l l2008bl ) . They calculated 
the evolution of the protostellar outflow only for a short duration of ^ 1 — lOyr at most, 
whereas the main accretion phase lasts for at least ^ 10^ — 10'^ yr. It is difficult to calculate 
the evolution of the protostellar outflow for a long duration by resolving the protostar itself 
because the timescales (time steps) in the regions around the protostar and the molecular 
cloud are quite different. 
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On the other hand, without resolving the driving source of the protostellar outflow, some 
studies focused on the long-term evolution of the protostellar outflow propagating into the 
interstellar m edium. In these studies, the outflow is input artificially in the computational 



domain (e.g.. lArce et al.l 120071 ). In addition, the circumstellar disk, the configuration, and 



the strength of the magnetic field are also adopted arbitrarily. Even with this type of 
calculation, we cannot estimate star formation efficiency because the mass ejection rate is 
artificially assumed in such studies. 



At the expense of spatial resolution around the protostar, iMachida et al.l (l2009al ) cal- 
culated the evolution of the collapsing cloud from the prestellar stage until the end of the 
main accretion phase. In their study, the driving source of the outflow (i.e., the first core and 
the circumstellar disk) was well resolved spatially, whereas the protostar was not resolved 
and was replaced by sink cells. They pointed out that the protostellar outflow reduces star 
formation efficiency to ~ 20 — 60%. However, they investigated only the evolution of a 
low-mass cloud core (M ~ 0.22 M©) and did not invest i gate the evolution of a cloud with 
a typical mass scale of ~ 1 Mq (e.g^ 



Motte et al 



1998 



Onishi et al.ll2002f ). In this study. 



using conditions similar to those in iMachida et al.l (j2009al ) , we calculate the evolution of a 
cloud with various initial masses and investigate the evolution of the outflow and the impact 
of the protostellar outflow on star formation efficiency. This paper is structured as follows. 
The framework of our models and the numerical method are described in §2. The numerical 
results are presented in §3. The mass ejection rate and star formation efficiency are discussed 
in 54. We summarize our results in §5. 



2. Model and Numerical Method 

2.1. Basic Equations 

To investigate star formation efficiency and long-term evolution of the protostellar out- 
flow, we calculate the star formation process from the prestellar core stage until the end of 
the main accretion phase for clouds with various masses using three-dimensional resistive 
MHD equations, including self-gravity: 

^ + V-(pt;) = 0, (1) 

dv 1 
p— + p{v ■V)v = -VP - —B X {y x B) - pV(t), (2) 
at Att 

f) Ft 

= V X (v X B) + r]V^B, (3) 

ot 

= AnGp, (4) 
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where p, v, P, B, rj, and (f) denote the density, velocity, pressure, magnetic flux density, 
resistivity, an d gravitational potenti a l, resp ectively. To mimic the temperature evolution 
calculate d bv iMasunaga fc Inutsuka I (120001). we adopt the piece-wise polytropic equation of 



state (see 



Vorobvov fc Basulbood : iMachida et al.lboO?! ) as 



-14 —3 

g cm 



2/3' 



rir 



(5) 



10^°cm-3). Equation 



shows 
^. For 



where Cs^o = 190 ms~^ and pc = 3.84 x 10" 

that the gas behaves isothermally for n ^ 10^° cm~'^ and adiabatically for n ^ 10^° cm" 
the realistic evolution of the magnetic field d uring star formation, we adopt a resistivity {rj) 
of the fiducial value in iMachida et al.l (120071 ). in which Ohmic dissipation be comes effective 



for 1 0^^ cm ^ ^ 10^^ cm ^ (for details, see Eqs. [9] and [10] and Fig. 1 of IMachida et al. 
20o3). 



2.2. Initial Settings 

This study investigates star formation efficiency in a single cloud. Thus, we assume 
an isolated cloud core embedded in the interstellar medium. As the initial state, we take a 
spherical cloud with a critical Bonnor-Ebert (BE) density profile pbe, in which a uniform 
density is adopted outside the sphere {r > Rc, where -Rc is the critical BE radius) to mimic 
the interstellar medium. The gravitational force is ignored outside the BE sphere (r > R^) 
to avoid the inflow of gas from outside the isolated core (i.e., from the interstellar medium). 
In addition, to strictly avoid mass inflow from outside the core, we prohibit mass inflow at 
r = Rc. Note that we do not prohibit mass outflow at the boundary between the BE sphere 
and the interstellar medium; the gas escapes freely from the BE sphere by the protostellar 
outflow. Hereafter, we call the gravitationally bound gas cloud within r < R^ the host cloud. 
We confirmed that the total mass of the host cloud is well conserved during the calculation 
before the protostellar outflow reaches the cloud boundary. Note that the total mass inside 
the host cloud is not conserved after the protostellar outflow propagates into the interstellar 
space (r > _Rc), because the mass is ejected from the host cloud by the protostellar outflow. 

Since the critical BE sphere is in equilibrium, we increase the density by a factor of / to 
promote contraction, where / is the density enhancement factor that represents the stability 
of the initial cloud. The density profile of the initial cloud is described as 

nfrl = / ^^^^""^ + ^ for r < , . 

' \ 0.02 pMRc) (1 + S,) f for r > R,, ^ ' 
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where pBE{r) is the density distribution of the critical BE sphere, and Sp is the axisymmetric 
density perturbation. An initial cloud with larger / is more unstable against gravity. The 
cloud stability is generally represented by a parameter an (= Et/Ep), which i s the r atio of 



thermal (E^) to gravitational (Eg) energy. As shown in iMatsumoto fc Hanawal (120031 ) . when 
the BE density profile is adopted, the density enhancement factor is related to the parameter 
ao as 

0.84 

"0 = (7) 

We adopt a density enhancement factor as 1.68, which corresponds to ao = 0.5. The density 
contrast between the center of the cloud (r = 0) and the cloud boundary (r = Rc) is 
p(r = 0)/p(r = Rc) = 14. In addition, a uniform density of pamb = 0.02pc (i.e, 2% of the 
central cloud density) is adopted outside the sphere (r > Rc). 

To break the axial symmetry, we add a small amount of m = 2-mode non-axisymmetric 
density perturbation to the initial core. For the m = 2-mode, in equation ([6]), we choose 

6p = A4r/Rc)\os2(P, (8) 

where (=0.01) represents the amplitude of the perturbation. The radial dependence is 
chosen so that the density perturbation remains regular at the origin (r = 0) at one time- 
step after the initial stage. This perturbation develops into a non-axisymmetric perturbation 
in the circumstellar disk that contributes to angular momentum transfer. In addition, this 
m = 2 perturbation ensures that the center of gravity is always located at the origin. 

For a dimensional BE density profile, we adopt an isothermal temperature of T = 10 K 
and a central number density of Uc = nc,o, where nc,o is a parameter in the range no = 
3 X 10^ — 3 X 10^ / cm~^. Since the temperature of each initial cloud is fixed, the size (i.e., 
radius and mass) of the BE sphere is uniquely determined only by the parameter n^^o- Thus, 
each model is characterized only by the initial central density, n^o- With these parameters, 
the critical BE radius (or radius of the host cloud) is i?c = 87 — 8700 AU. The mass inside 
r < Rc for each model is Md = 0.015 — 1.5 Mq. The host cloud radius and mass for each 
model are listed in Table [TJ 

In each model, the cloud rotates rigidly around the z-axis in the r < Rc region and 
a uniform magnetic field parallel to the 2;-axis (or rotation axis) is adopted in the entire 
computational domain. The magnetic field strength and the rotation rate are scaled using 
the central density po = PBE^r = 0)/ as 

a = i?oV(4vrpo<o), (9) 
u = no/{^7rGpo)'/^. (10) 
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In all models, we adopt a = 0.1 and u = 0.1, whic h are the most appropr i ate parameters 



for dr iving; strong outflow in the collapsing cloud ( iMachida et al. I l2005bl : iMachida et al. 
2008bl ). Since the magnetic field and rotation are normalized by the central density and 



each model is characterized by only the central density, each cloud has a different magnetic 
field strengths and angular velocity. The magnetic field (Bq) and angular velocity (Qq) 
for each model are summarized in Table [H However, all models have the same ratios of 
rotational and magnetic energies to the gravitational energy, /3o (= -ErotZ-E'grav = 5 x 10~^) 
and 7o (= -EmagZ-^'grav = 4 X 10"^), where E^^t and -Emag are rotational and magnetic energies, 
respectively. In addition, all models have the same mass-to-fiux ratio M/$. There exists 
a critical value of M/$ below which a clo ud is supported against gravit y by the magnetic 
field. For a cloud with uniform density, iMouschovias &: Spitzerl ( 119761 ) derived a critical 
mass-to-fiux ratio 

where the constant C, = 0.48 f Tomisaka et al. IQSSal jbl). The mass-to-fiux ratio normalized 
by the critical value A is described as 



1/2 



:ii^ 



:i2) 



All models have A = 4, which is slightly larger than the typical value of observation. The 
observations indicate that molecular cloud c ores have the rn ass-to-fiux ratio in the range of 
0-8 ;^ A ;^ 7.2 with a median value of A ^ 2 ( jCrutcheil 119991 ). 



The model names and parameters are also listed in Table [H We believe that these energy 
ratios are adequate for comparing cloud evolution among models with different masses. We 
can match the magnetic field strength (Bq) and rotation rate (Qq) for clouds with different 
masses. However, since doing so changes the evolution of the cloud and the protostellar 
outflow, it is difficult to compare the star formation efficiency of the models. For example, 
when the initial magnetic field strength is fixed among the models, no outfiow may appear in 
a less massive cloud because it has a small ratio of magnetic energy to gravitational energy. 

In this paper, to compare cloud evolution among models, we usually use the freefall 
timescale at the center, r = (tff,c)) and boundary, r = Rc (%,b), of the initial cloud as the 
unit of time. Since the density contrast between the center of the cloud and its boundary is 
14, the freefall timescale of the cloud boundary is about 3.7 times longer than that at the 
center of the cloud (tfr^b = 3.7tff,c)- The freefall timescale at the center of the initial cloud 
(^ff,c) for each model is listed in Table [TJ 
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2.3. Sink Cell and Numerical Method 



To realize the long-term calculation of star formation, we adopt a sink at the center of 
the cloud. We start the calculation without a sink and calculate the cloud evolution for the 
prestellar gas collapse phase without a sink. Later, we identify protostar formation in the 
collapsing cloud core when the number density exceeds n > rit^r at the cloud center, where 
nthr is the threshold density. After protostar formation, we calculate the cloud evolution 
with the sink. 

To model the protostar, we adopt a fixed sink with a radius of 1 AU composed of sink 
cells only around the center of the computational domain. Since we add only m = 2 density 
perturbation, as described in §2.21 the protostar (or center of gravity) does not move and 
remains at the center of the computational domain during the calculation. In the region 
r < Tsink = 1 AU, the gas having a number density of > rithv = 10^^ cm~'^ is removed from 
the computational d omain and added to t he protostar as a gravitating mass in each time 
step (for details, see iMachida et al.l l2009al ). Thus, for each time step, the accretion mass 
onto the protostar is calculated as 



[p{hj, k) -pthr] dV. 



(13) 



^■<''smk 



In addition, inside the sink, the magnetic flux is removed by Ohmic dissipation because 
this region has a na agnetic Reynolds number Re exceeding unity Re > 1 (for details, see 
Machida et aPboOTh . 



see 



For calculation on a la rge spatial scale, the nested grid method is adopted (for details. 



Machida et al. ll2005al Jbl). Each level of a rectangular grid has the same number of cells 
(64 X 64 X 32). The calculation is first performed with five grid levels (/ = 1 — 5). In all 
models, the fifth level of the grid (/ = 5) has a box size of = 2 Rc and just covers the 
entire region of the isolated BE sphere. The first level of the grid has a box size of Li = 2^ R^. 
and is filled with low-density interstellar medium outside r > R^. Thus, we can calculate 
the propagation of the protostellar outflow in the region of < 2^ R^. The protostellar outflow 
never reaches the computational boundary by the end of the calculation. In addition, we have 
checked that the Alfven waves generated at the center of the cloud (or the computational 
boundary) never reaches th e computational bound ary (or the center of the cloud) during the 
calculation (for details, see iMachida et al.ll2010bl ). After the calculation starts, a new finer 
grid is generated before the Jeans condition is violated (ITruelove et al. 1119971 ). Although 
the maximum grid level differs among the models, each model has a spatial resolution of 
< 0.3 AU in the finest grid. 
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3. Results 



In this section, we present the evolution of the collapsing cloud and the outflow for 
typical models in §3.1.11 (model N08, less massive cloud case) and §3.1.21 (model N06, massive 
cloud case). Then, we describe the mass accretion rate ( §3.2p and masses of the protostar 
and circumstellar disk ( §3.3p . and the mass ejected by the protostellar outflow ( §3.4p . The 
properties of the protostellar outflow are described in §3.51 



3.1. Typical Model 



3.1.1. Model N 08 



Figure [T] plots the cloud evolution for model N08 from the initial state until the end of 
the main accretion phase; only the z > region is shown. Note that although the structure 
only around the initial cloud scale (~ i?c) is plotted in Figure [H the computational domain 
has a size of 2^i?c- As seen in Figure [T]a, we adopted a spherical cloud with the BE density 
profile and the radius of Rc as the initial state. As described in §2.2[ since we ignored 
the gravity outside the BE sphere (r > i?c), only the gas in the r < region (inside 
the white dotted line in Fig. [1]) can collapse to fall onto the center of the cloud. For this 
model, the number density exceeds n > rithr and a protostar forms at t = 3.51tfj^c, where tg^c 
(= 1.4 X 10^ yr) is the freefall timescale of the initial cloud at the center. Figure [T]& shows 
the density distribution just after protostar formation. The red contours of n = 10 cm~^ in 
Figure [T]a and b indicate that the cloud gradually contracts toward the center with time. 

The blue contour in the figure corresponds to the boundary of the protostellar outflow, 
inside which the gas moves outward against the center of the cloud. Thus, this contour shows 
the shape of the protostellar outflow. The outflow is driven by the circum stellar disk that 



originates from the first core, which appears before protostar formation (ITomisakal 12002 



Machida et al. I l2005bl : iHennebelle fc Fromang] l2008al : iBiirzle et al.1 120111). After protostar 



formation, the first core becomes the circumstellar disk ( Bate 



1998 



Inutsuka et al. 2010 



Machida et al.ll2010al ). The outflow is then driven by th e circumstellar disk af t er a s mooth 
transition from the first core to the circumstellar disk (IMachida et al.l l2009ai 1201 if ) . The 
protostellar outflow remains in the host cloud (r < Rc) for t ^ 7000 yr (Fig. [T]c); it reaches 
the boundary of the host cloud at t ~ 7200 yr (Fig. [T](i) and penetrates the host cloud 
boundary and propagates into the interstellar space (r > R^) for t ^ 7200 yr (Fig. [T]e-/i). 
The outflow driving halts inside the host cloud at t ~ 17000 yr (^ 12tfj c)- The infalling 
envelope is depleted and the mass accretion is almost over by this epoch ( §3.30 . Thus, the 
circumstellar disk cannot drive the outflow at this epoch because the outflow is powered by 
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gas accretion onto the circumstellar disk. Since the freefall timescale of the cloud boundary 
is tfT,b = 5300 yr, a large fraction of the gas inside the host cloud has already fallen onto 
either the circumstellar disk or the protostar by this epoch. Although the outflow at the 
host cloud scale completely disappears in t ~ 3 x 10"^ yr (=20tff^c = 5.4 tg b), the density 
cavity formed by the protostellar outflow remains around the host cloud as seen in Figure [T]j. 

Figure [1] also shows that the infalling envelope in the host cloud depletes with time. 
Part of the gas inside the initial host cloud is ejected by the protostellar outflow, and the 
remainder falls onto either the circumstellar disk or the protostar. In Figures dFi and i, 
we can see only a disk-like structure with the size of ~200AU, because the infalling gas is 
depleted by this epoch. At these epochs, the density of the infalling envelope [n ^ 10^ cm~^) 
is less than 1/100 the initial cloud density {n ^ lO^cm"^). At the end of the calculation 
[t ^ 20tff c), the mass ratio of the infalling envelope to the total mass of the initial cloud is 
only < 4%. Thus, the main accretion phase has already ended by this epoch. 

The evolution of the outflow configuration is shown in Figure [2l Each panel corresponds 
to the same epoch as in Figure [T](i, e, /, and h. Figure [2] shows that the protostellar outflow 
propagates into the interstellar space while maintaining good collimation. However, Fig- 
ure [2]a and b indicate that, early in the evolution, the outflow also extends in a transverse 
direction and increases in width. After its width becomes comparable to the host cloud 
scale, the outflow extends only in the vertical direction while maintaining its width. The 
magnetic field lines that guide the outflow are anchored by the host cloud (or the gravita- 
tionally bound sphere). Thus, the final width of the outflow is comparable to that of the 
host cloud. Therefore, the outflow collimation improves with time after the outflow width 
becomes comparable to the host cloud scale, whereas the collimation is not good when its 
width is smaller than the host cloud scale or the outflow remains within the host cloud. 



3.1.2. Model N 06 

Figures [3] and H] show the density and velocity distributions for model N06 at t = 
1.315 X 10^ yr (= 9.4tff^c=2.5tff,b)- For this model, the protostar forms at t = 4.502 x 10'^ yr 
after the cloud begins to collapse. Thus, the figures show the structure 8.647 x 10"^ yr after 
the protostar formation. By this epoch, the outflow penetrated the host cloud that has 
a radius of 4800 AU and reaches ~ 2 x 10^ AU from the center of the cloud, as seen in 
Figure [3] left panel. The figure also shows the bow shocks caused by the protostellar outflow 
at 2 ~ ±2 X 10^ AU. The upper right panel in Figure [3] plots the structure around the host 
cloud and shows that the gas flows out from the host cloud in the vertical direction by the 
protostellar outflow. In addition, inside the white dotted line that denotes the size of the 
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initial cloud, the gas density is considerably lower than the ambient medium. This is because 
a part of the gas falls onto the central region to form the protostar and circumstellar disk, 
while the remainder is ejected from the host cloud by the outflow. The lower right panel is 
16 times magnified view of the central region of the upper right panel. This panel shows that 
the outflow is strongly driven by the disk-like structure around the protostar. On this scale, 
the gas accretes onto the protostar through the circumstellar disk, and a part of the accreting 
mater is ejected by the outflow. Figure H] shows the structure around the circumstellar disk 
on y = (left panel) and z = (right panel) planes. These panels show that the rotating 
disk exits around the protostar and drives the outflow. Thus, it is expected that th e outflow 
is mainly driven by the magnetocentrifugal mechanism (IBlandford fc Paynelll982l ). 

Figures [3] and H] indicate that large scale outflow with a size of > lO'* AU (Fig. |3] left 
panel) originates from the disk wind driven by the circumstellar or rapidly rotating disk with 
a size of ~ 1 — 100 AU. The disk wind propagates into the infalling envelope, and thus it 
sweeps and collects a larger fraction of the infalling matter. Finally, a large fraction of the 
swept material is expelled from the host cloud. Thus, the outflow significantly affects the 
final protostellar mass or star formation efficiency. 



3.2. Mass Accretion Rate onto Protostar 

Figure [5] shows the mass accretion rate onto the protostar and the protostellar mass 
against the time after protostar formation t. Here, we describe the time after the protostar 
formation as t, which is defined as 

i=t-to, (14) 

where t is the elapsed time after the cloud begins to collapse (or the calculation starts), and 
to is the protostar formation epoch and is listed in Table [H In each model, the gas density 
exceeds n > nthr, and the protostar forms ~ 3 — 5tff,c after the cloud begins to collapse 
(Tabled]). Note that model N39 shows no continuous collapse and no protostar appears; this 
is because the initial cloud density for model N39 is too high to induce continuous collapse. 
In our models, since the gas becomes adiabatic at n = ric ^ 10^° cm~^ ( §2.21) . the initial cloud 
density {ric = 3 x 10^ / cm^^) for model N39 is very close to this critical value ric- As a result, 
model N39 shows repeated contraction and expansion around the initial configuration, not 
continuous collapse. Below, we describe only the models showing protostar formation (N35, 
N06, N36, N07, N37, N08, N38, N09) and do not mention model N39 again. 

As described in §2.31 we removed the gas having a number density of n > rithr in 
the region r < 1 AU from the computational domain. We regarded the removed gas as 
the accreted mass onto the protostar and estimated the mass accretion rate in each time 
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step. Figure [5] shows that, in each model, the mass accretion rate is as high as M ~ 
10~^ — 1O~^M0 yr~^ just after protostar formation {t ^ Oyr). In theoretical analyses, the 
mass accretion rate is defined as M = f c^ /G, where / is a numerical factor (e.g., / = 



0.975 in iShul 119771 . / = 46.9 in iHunterl 119771 ). Since gas clouds have temperatures of T = 
10 K (cs = 0.2kms~^), the accretion rate in the main accretion phase is M = (2 — 90) x 
lO~^M0yr^^. Thus, the accretion rate derived in our calculations corresponds well to the 
theoretical expectation. 

Then, in each model, the mass accretion rate gradually decreases with time, to M ~ 
10"^ - 1O-^M0 yr"^ at t ~ ts,c and M ~ 10"^ - lO-^M© yr'^ at t ~ ts,h- Thus, gas accretion 
almost halts and the protostar rarely increases in mass at t ^ tg^b- As shown in Figure [5], all 
models show a qualitatively and quantitatively similar mass accretion rate trend when the 
mass accretion rate onto the protostar is normalized by the freefall timescale of the initial 
cloud. However, since the real (or dimensional) time of the freefall timescale depends on the 
initial cloud density (or the initial cloud mass), the duration of the main accretion phase 
is different in models with different cloud masses. A cloud with a larger mass (or lower 
density) has a longer duration of the main accretion phase. For example, model N35 with 
Mci = 1.5 Mq has the mass accretion of M > lO~^M0yr~^ for ~ 5 x 10^ yr, while model 
N09 with Mci = 0.03 Mq has that only for ~ 10'^ yr. Figure [5] also shows that the protostar 
formed in a massive cloud is more massive than that formed in a less massive cloud. This 
is because the massive cloud has a longer gas accretion phase, so the protostar has enough 
time to acquire sufficient mass and evolves into a relatively massive star. 

In Figure [5], the mass accretion rate oscillates in models having an initially massive 
cloud. This is caused by a non-axisymmetric structure appearing in the circumstellar disk. 



(Bate 


1998; 


Walch et al. 


2009a: 



axisymmetric or spiral structure develops owing to gravitational instability. Such structure 
effectively transfers angular momentum outward and intermittently promotes mass accretion 
onto the protostar. As a result, these models sh ow time variability in the mass accretion 



rate flVorobvov fc Basull2006l : Machida et al.ll201lh . In Figure[5l models N35, N06, N36, N07 



and N37 show time variability in the mass accretion rate. 



3.3. Mass Evolution of Protostar and Circumstellar Disk 

The mass of the protostar, circumstellar disk, protostellar outflow, and infalling envelope 
are plotted against time after protostar formation t in Figure [61 In addition, these masses at 
iff after protostar formation are listed in Table [H The circumstellar disk mass is estimated 
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according to the formula in lMachida et al.l ( l2010al ). To estimate the outflowing gas Mout, we 
integrated the gas with velocity Vr > for the entire computational domain and subtracted 
the gas swept by the protostellar outflow outside the host cloud r > R^. from the integrated 
mass. For later convenience, we divide the outflowing mass into two parts, 



(15) 



where Mgj is the mass ejected from the host cloud, and Mout.Rc is the outflowing mass having 
Vr > Cs inside the host cloud (r < Rc). To calculate the mass of the infalling envelope, we 
calculated the total mass Mtot inside the host cloud (r < Rc). Then, we subtracted the disk 
mass Mdisk and the outflowing mass Mout,Rc in the r < Rc region from the total mass. 



Mtot(r < Rc) - Mdisk - Mont,RAr < Rc 



(16) 



In our models, the sum of the protostellar mass, circumstellar mass, and infalling envelope 
mass is not conserved inside the host cloud (i.e., r < Rc) because part of the gas is ejected 
from the host cloud by the protostellar outflow. However, we confirmed that the sum of 
the total mass and the protostellar mass is well conserved before the outflow reaches the 
boundary of the host cloud, as described in §2.21 

Figure [6] shows that, in each model, the mass of the infalling envelope decreases to 
Menv/Mci ^ 0.1 at t ~ tfj^b (Table [2]), where is the initial cloud mass. Thus, the mass of 
the infalling envelope is depleted and the main accretion phase is almost over by this epoch. 
Note that, in Figure El we subtracted the time (to) until the protostar forms from the time 
(t) after the cloud begins to collapse. Thus, at the epoch ts^h, indicated by the arrow in 
Figure [6l a longer time than tg^b has passed since the cloud begins to collapse (Table [1]). 
Therefore, it is reasonable that almost all the gas has already fallen onto the center of the 
cloud at this epo ch i ~ tff^hi because th e gas falls onto the cloud center in several times 
freefall timescale ( iMachida et al. I l2005al ) and the freefall timescale at the cloud boundary 
(tff^b) has already passed by this epoch. 

Figure E] also indicates that the circumstellar disk mass dominates the protostellar mass 
just after protostar forma tion in each model. The circumstellar has been reported to orig- 
inate from the first core ( llnutsuka et al.l 120101 ). which is about 10- 1 00 times more massive 



than t he protostar at the pro tostar formation epoch (ILarsonl Il969l : iMasunaga fc Inutsuka 



2000l ). iMachida et al.l (j2010al ) showed that the circumstellar disk is inevitably more massive 
than the protostar in the early main accretion phase because the first core evolves directly 
into the circumstellar disk after protostar formation. In the later main accretion phase, 
however, in some models, the protostellar mass dominates the circumstellar disk mass. In 
the initially massive clouds (models N35, N06, and N36), the protostellar mass continues to 
increase, whereas the circumstellar disk mass gradually decreases until the end of the main 
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accretion phase. In addition, in models N35 and N06, the protostar becomes more massive 
than the circumstellar disk for t ^ c- The massive circumstellar disk becomes gravitation- 
ally unstable and tends to exhibit a non-axisymmetric structure that can transfer angular 
momentum outward; the gas in the circumstellar disk effectively falls onto the protostar. 
Therefore, the accretion rate onto the protostar from the circumstellar disk dominates the 
accretion rate onto the circumstellar disk from the inf ailing envelope, and the mass of the 
disk begins to decrease. 

Figure [7] shows the density distribution on the equatorial plane around the center of 
the cloud at t ~ 7ts,c for models N35, N06, N08 and N09. The figure shows that a rotating 
disk with a size of ~ 10 — 100 AU forms around the center of the cloud by this epoch. 
As seen in Figure [7]a, two clumps appear in model N35, in which fragmentation occurs in 
the circumstellar disk about ~ 5 x 10^ yr after protostar formation. In models N35 and 
NOG, the circumstellar disk is rather massive in the early main accretion phase. As shown 
in Figure [6l the circumstellar disk has a mass of Mdisk ~ 0.4 — 1.0 Mq in model N35 and 
~ 0.4 Mq in model NOG in the early main accretion phase for t <^ tg^c- Such a massive 
disk becomes gravitationally unstable and tends to develop a non-axisymmetric structure 
and show subsequent fragmentation, as seen in Figure [7]a. Although model NOG shows 
no fragmentation until the end of the calculation, non-axisymmetric (or spiral) structure 
develops in the circumstellar disk, as seen in Figure [7]?). The circumstellar disk mass begins 
to decrease after the non-axisymmetric structure develops because such structure transfers 
angular momentum outward, and the gas in the circumstellar disk effectively falls onto the 
central protostar. Finally, the protostellar mass dominates the circumstellar disk mass in the 
later main accretion phase for models N35 and NOG (Fig. [6]). For model N3G, the circumstellar 
disk mass gradually decreases in the main accretion phase (Fig. [6]), whereas it is greater than 
the protostellar mass at t = tfj.b- The efficiency of angular momentum transfer for model 
N3G is considered to be lower than those for models N35 and NOG because only a weak spiral 
structure appears in this model. 

On the other hand, for models N07, N37, N08 and N38, the circumstellar disk mass does 
not decrease greatly in the main accretion phase and is slightly greater than or comparable to 
the protostellar mass at the end of the main accretion phase. Although the circumstellar disk 
is more massive than the protostar during the main accretion phase, no non-axisymmetric 
structure develops in these models. This is because the size of the circumste llar disk is not 



sufficiently larger than the Jeans length, and the disk is stable against gravity (IMachida et al. 



2010al ). As shown in Figure [71 model NOB has an almost axisymmetric structure. Thus, 
angular momentum transfer due to non-axisymmetric structure is not so effective, and a 
massive disk remains until the end of the main accretion phase. 
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As seen in Figure [6|, for model N09, the protostar becomes more massive than the 
circumstellar disk just after the protostar formation, and the protostellar mass dominates 
the circumstellar disk mass by the end of the calculation. Since model N09 has no sufficient 
cloud mass at the initial state, the circumstellar disk cannot increase its mass significantly by 
gas accretion. For this model, the protostar and the circumstellar disk just after the protostar 
formation has a mass of Mps ~ 0.01 and Mdisk ~ 0.01 M0, respectively. Thus, the sum 
of protostellar and circumstellar disk masses (~ 0.02 M©) is comparable to the initial cloud 
mass (Mci = 0.03). Therefore, mass accretion onto the circumstellar disk from the infalling 
envelope rapidly declines immediately after the protostar formation. On the other hand, the 
protostar gradually increases its mass by the mass accretion from the circumstellar disk. 

Figure E] indicates that, in each model, the protostellar mass rapidly increases for t < tg_c, 
and slightly increases for ^ ^ b- For i ^ tg b, the protostellar mass rarely increases 
because the infalling envelope is almost depleted by this epoch. To compare the protostellar 
mass evolution among models, the time averaged mass accretion rate for i < tg^h is estimated 
as Mave = Mps(tfr,b)/^ff,b, where Mps(tfr,b) is the protostellar mass at t = tg^b- For all models, 
the averaged mass accretion rates are in the range of 4.0 x 10~^ < Mave/(MQ yr~^) < 
8.5 X 10~^. The difference of the mass accretion rate among models is only a factor of 
about 2. Thus, a protostar has a similar accretion history when the protostellar evolution is 
normalized by the freefall timescale of the initial cloud. In reality, however, initially massive 
cloud has a longer freefall timescale and longer duration of the main accretion phase. Thus, 
the protostar formed in massive cloud has a greater mass at the end of the main accretion 
phase. 



3.4. Mass Ejected by Protostellar Outflow 

In the main accretion phase, the mass of the protostellar outflow dominates, or is com- 
parable to, both the protostellar and circumstellar disk masses for models N35, N06, N36, 
N07, and N37, as shown in Figure Ei This indicates that the protostellar outflow greatly af- 
fects star formation efficiency because it ejects a large fraction of the mass of the host cloud. 
In other words, the protostellar outflow controls t he protostellar mas s or st ar formation effi- 



ciency in the star formation process as pointed by lMatzner fc McKed (120001 ). Figure [8] shows 
the density and velocity distribution on the y = cutting plane for models N35, N06, NOB, 
and N09 at t ~ Ttg c ~ 2tg^b- Note that the models and epochs in Figure M are the same 
as those in Figure [71 but the spatial scales differ. In each panel, the outflowing region is 
denoted by a white contour, inside of which radial velocity of the gas is supersonic {vr > Cg). 

Figure |H] shows that the gas is strongly ejected from the host cloud by the protostellar 
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outflow in any model. In addition, the outflow width (i.e., the horizontal length in the 
direction of travel) is comparable to the radius of the host cloud (Re), and thus, the outflow 
has a wide opening angle. At this epoch, on the cloud scale of r = Rc, the outflow has 
opening angles of 42° (N35), 45° (N06), 37° (N08), and 26° (N09). Note that the opening 
angle continues to increase until it becomes comparable to the cloud radius. As a result, 
the protostellar outflow with a wide opening angle sweeps up a large fraction of the infalling 
material and ejects it into the interstellar space. In addition to this swept material, a part 
of the mass in the circumstellar disk is ejected directly by the outflow. 

Figure |9] shows the ratio of the outflowing to infalling masses for models N35, N36, 
N37, and N38. To investigate outflow efficiency and mass ejection rate, we calculated the 
outflowing/inf ailing mass rate on the / = /max — 1 and /max — 3 grid surfaces. In each model, 
the / = /max — 1 gi'id covers the entire circumstellar disk, and / = /max — 3 grid has a size 
of ~ 8 — 10 times the disk radius. Since the protostellar outflow is originally driven by 
the circumstellar disk, the mass ejected from the / = /max — 1 grid almost corresponds to 
the outflowing mass directly driven by the circumstellar disk. On the other hand, the mass 
ejected from the / = /max — 3 is the sum of the outflowing mass directly driven by the 
circumstellar disk and the swept mass b y the outfl o w tha t propagates into the dense (or 



massive) infalling envelope. According to iTomisakal (|2002[ ). we estimated the outflow Mout 
and inflow Min masses as 



-^out = / pma.x[v ■ n, 0] ds, (17) 

' boundary of 1 



and 



Min = / pmsix[v ■ —n,0]ds, (18) 

J boundary of 1 

respectively, where n is the unit vector outwardly normal to the surface of the / = /max ~ 1 
or /max ~ 3 grid. Figure M shows that the ratio of the outflowing to inflowing mass around 
the circumstellar disk (/ = /max — 1, thin line) has a peak of Mont/Mi^ 1 at t ^ ts^c- 
Thus, the mass ejection rate is comparable to the mass infalling rate at this epoch. Note 
that since the circumstellar disk and the outflow driving region are embedded in the grid of 
/ = /max ~ 1; the mass swept up by the outflow may be slightly included in the estimation 
of the outflowing rate. Then, in all the models, the ratio in the / = /max ~ 1 grid gradually 
decreases with time. The ratio decreases to Mout /Af in ~ 0.1 — 0.3 at t ~ tg^b and reaches 
Mout/Min < 0.01 at t ~ 10tfr,c. 

Figure also shows that the ratio of Mout/^in in the hmax — 3 grid exceeds that in the 
/max ~ 1 grid for a short period after the protostar appears. Thus, the outflowing mass (rate) 
increases with distance from the driving source. This indicates that a wide-opening-angle 
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outflow sweeps up the infalling gas and ejects it toward the center of the cloud. Especially, for 
t ^ tfj c, the rate in the / = /max — 3 grid is about 10 times larger than that in the / = /max — 3 
grid. This indicates that the protostellar outflow collects the matter in the infalling envelope 
10 times more massive than that directly blown away from the circumstellar disk. Note that 
since the boundary of the host cloud is located more far away from / = /max — 3 grid, the 
outflow sweeps more matter in the infalling envelope to be ejected from the host cloud. 

The thick line in Figure [H] indicates that the outflow on a larger scale gradually decreases 
for t ;^ 1 — Stff^c and disappears in t ^ lOtg^c- Thus, the lifetime of the outflow is about 
10 times the freefall timescale of the center of the host cloud, or 3-5 times of the freefall 
timescale of the outer edge of the host cloud. The protostar forms ~ 3 — 5 tg^c after the cloud 
begins to collapse (Tabled]). Thus, roughly speaking, the protostellar outflow continues to 
be driven by the circumstellar disk until the freefall timescale of the cloud boundary passes 
after the protostar formation. This is natural that the almost all the gas inside the cloud fall 
onto the circumstellar disk in b (Figs. [5] and [6]), and the protostellar outflow is powered 
by the mass accretion onto the circumstellar disk. 



3.5. Evolution of Protostellar Outflow and its Collimation 

Figure dU] shows the outflow length (upper panel) and width (lower panel) against the 
time normalized by the freefall timescale at the center of the cloud for all models. Note 
that the length and width in Figure [10] right panels are normalized by the initial each cloud 
radius Re- We defined the outflow as the gas having the (positive) radial velocity larger 
than the sound speed [vr > Cg). The outflow expands with time and reaches ~ 500 — 10^ AU 
at ~ lOtff^c (Fig. [TU] upper left panel). The flnal size of the outflow depends on the size 
(or mass) of the initial cloud. The outflow extends up to about 10 times the initial cloud 
radius except for model N09 (Fig. [10] upper right panel). The outflow in a massive cloud 
has a longer lifetime to reach a greater distance from the protostar because a massive cloud 
has a longer freefall timescale and the outflow continues to be driven on ~ 1 — 10 times the 
freefall timescale. The lower panels of Figure [10] indicates that the outflow also expands in 
the horizontal direction and becomes comparable in size to the host cloud radius Rc. The 
outflow propagates into the infalling envelope along the magnetic fleld lines. Although the 
magnetic fleld lines that drive the outflow are strongly bundled around the circumstellar 
disk, they spread with the distance from the center of the cloud up to ~ Rc because they 
are connected to the cloud scale lines. Therefore, the outflow width also spreads in the 
horizontal direction and has an opening angle of ~ 45° on the host cloud scale. 

Figure [TT] shows the evolution of the collimation factor of the outflow, which is deflned 
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as 

/col = (19) 

Wout 

where rout and Wout are the length and width of the outflow, respectively. The figure indicates 
that the coUimation factors remain /cooi ;^ 10 for t < 5 — 7tfj.c- Then, they begin to increase 
for t ^ 5— Ttff^c- At the end of the calculation, the coUimation factors reach /cod ~ 10—30. As 
seen in Figure [2l although the outflow has a wide-opening-angle in the early main accretion 
phase, it is a very well collimated in the later main accretion phase. In addition, the width of 
the protostellar outflow reflects the size of the host cloud. Thus, we can acquire information 
on the prestellar cloud core from the size and width of the protostellar outflow. 



4. Discussion 

4.1. Outflow Momentum and Comparison with Observations 

To investigate the evolution of the outflow strength, we estimated the outflow momen- 
tum, which is defined as 

MKut = J pvontdv, (20) 

where fout is the outflow velocity. Figure [12] upper panel shows that the evolution of the 
outflow momentum against the elapsed time after the outflow appears tout? which is defined 

as 

^out = t — tout.O) (21) 

where tout,o is the time at the moment of the outflow appearing. This panel indicates that an 
initially massive cloud has a larger outflow momentum. The outflow momentum increases for 
t ^ ifT,b after the outflow appears, whereas it gradually weakens for t ^ b- Since the main 
accretion phase is almost over at t ~ tg^b? the outflow momentum has a peak at t ~ tg ]^. 
At its peak, the outflow momentum is 0.003 — O.6M0kms~^, depending on the initial cloud 
mass. The outflow originating from an initially massive cloud has a larger momentum and 
reaches further, while that from an initially less massive cloud has a weak momentum and 
disappears in a short duration. 



Arce et al.l (120071 ) observed many protostellar outflows in the Perseus molecular cloud 
complex and statistically investigated them in detail. They showed that a large fraction 
of protostell a r out flows have momenta in the range of 0.05 ^ MV^ut/( ^©kms"^) ^ 1. 



Curtis et al.l (l2010l ) also investigated outflows in Perseus molecular cloud and showed that 
outflow momenta are distributed around MV^ut ~ O.lMokms"^. Figure [T2l upper panel 
shows that models with typical cloud mass of Md > 0.1 M0 (models N35, N06, N36, N07 
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and N37) have the outflow momentum of MV^ut ~ 0.01 — O.6M0kms~^ for tout ^ 10^ yr. 
Thus, the observation of outflow momentum well agrees with simulation results. Note that 
the observation show a snapshot of many protostellar outflows, while Figure [T2] upper panel 
shows the momentum evolution of individual outflow in each cloud. 



Figure [12] lower panel shows the momentum flux of the outflow, which is defined as 

^ MVout 



tout 

1-7 



(22) 



The figure shows that the protostellar outflow has a momentum flux in the range of 10^^ 
F/( M0kms~^/yr) ^ 10^'^. The momentum flux gradually decreases with time. Roughly 
speaking, the outflow has a momentum flux of F ^ 10"^ Mq kms~^/yr in the e arly main ac- 



cretio n phase and F ~ 10 ^ Mq km s ^/yr in the late stage of the star formation. iBontemps et al 



(]l996l ) showed that, with 45 observed outflow samples, the outflow momentum flux is typ- 



ically F ~ 10 ^ Mp km s ^ at the ear ly Class stage and F ~ 2 x 10 ^ Mokms ^ at the 



late Class I stage. ICurtis et al.l (|2010[ ) also shows the similar trend of the momentum flux. 



Thus, the momentum flux derived from our simulation well agrees with observations. 

In the calculations, we adopted the minimum spatial scale of ~ 0.3 AU as described 
in §2.21 Thus, we could not resolve protostar and its neighborhood. Therefore, no high 
velocity jet that is driven near the protostar appears. However, outflow momenta derived in 
our simulation well agree with observations. This indicates that observed outflow momenta 
can be explained only by the outflow driven by the circumstellar disk without both high 
velocity jet and mass entrained by the jet. In addition, since the outflow driven by the 
circumstellar disk has a wide-opening-angle, it greatly contributes to determine the star 
formation efficiency. 



4.2. Final Mass of Protostar, Disk and Protostellar Outflow 

The masses of the protostar, circumstellar disk, and outflow at t ~ to + %,b are listed 
in Table |21 The table also lists the mass fraction of the protostar (cps), protostar plus 
circumstellar disk (edisk), outflow (eout) and infalling envelope (ecnv) to that of the initial 
cloud at t ~ to + ^ff,b- As shown in Figure |5l since the mass accretion rate at this epoch 
(t ~ to + tfr_b) decreases to Mps <^ 1O~^M0 yr~^, the protostar cannot acquire additional mass 
from the infalling envelope after this epoch. Figure E] and Table [2] indicate that more than 
90% of the initial cloud mass is already depleted at this epoch. Part of the cloud mass is 
ejected from the host cloud by the protostellar outflow; the reminder falls onto either the 
circumstellar disk or the protostar. In summary, the main accretion phase is almost over by 
this epoch. 
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The upper panel of Figure [13] shows the masses of the protostar, outflow, and circum- 
stellar disk at the end of the main accretion phase (t ~ to + ^fr,b) against the initial cloud 
mass. The protostar and outflow increase in mass as the initial cloud mass increases. On 
the other hand, the mass of the circumstellar disk increases with the initial cloud mass when 
Mci < 0.5 Mq, whereas it saturates when Mdisk ~ 0.1 for Md > 0.5 Mq. As seen in 
Figure [TJ a massive protostar formed in a massive host cloud has a massive, gravitationally 
unstable circumstellar disk. Such a disk shows spiral structure and subsequent fragmen- 
tation. The spiral structure or fragments effectively transfer angular momentum outward 
and promote mass accretion from the circumstellar disk onto the protostar. Therefore, the 
massive circumstellar disk is self-regulated: when it is sufficiently massive, mass accretion 
onto the protostar is amplified, and the disk mass begins to decrease. 

The upper panel of Figure [12] also shows that the mass of the protostellar outfiow is 
larger than or comparable to the protostellar mass for Md > 0.08 M©. This indicates that, 
in a massive host cloud, the protostellar outfiow sweeps up a large fraction of the infalling 
mass and ejects it into the interstellar space. Thus, the protostellar outfiow significantly 
affects star formation efficiency in such a cloud. On the other hand, in a less massive cloud, 
the mass of the protostellar outflow is smaller than the protostellar mass. As described in 
^ the flrst core is formed before the protostar formation. At its formation, the flrst core 
has a mass of ~ 0.1 — 0.01 which is comparable to the Jeans mass at this epoch. The 
flrst core or the circumstellar disk drives a wide-opening-angle outflow. However, when the 
initial cloud mass is comparable to the mass of the flrst core, the infalling envelope does not 
contain enough mass to be swept up by the protostellar outflow. As a result, the outflow 
cannot accumulate enough mass and only a small fraction of the host cloud mass is ejected. 
Thus, higher star formation efficiency is realized in a less massive cloud without significant 
mass ejection from the host cloud. 

Figure [T3] lower panel shows the mass ratio of protostar, circumstellar disk, and proto- 
stellar outfiow at t = to + ^fr,b- As shown by the blue line, the circumstellar disk mass is 
comparable to the protostellar mass for Md < 0.26 M©; in this case, the circumstellar disk 
has an almost axisymmetric structure. On the other hand, for M^i > 0.26 M©, the mass 
ratio of the circumstellar disk to the protostar decreases rapidly. The mass ratio of disk-to- 
protostar for model N09 is Mjisk/Mps = 0.71, and that for model 35 is M^^sk/Mps = 0.22. 
The difference is the result of the non-axisymmetric structure and fragmentation appearing 
in the circumstellar disk. Such structure greatly promotes mass accretion from the circum- 
stellar disk onto the protostar and decreases the circumstellar disk mass. The cloud with 
Mci ~ 0.1 — 0.5 Mq has a marginally gravitationally stable disk and has a maximum mass 
ratio of the circumstellar disk to the protostar. 
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The black line in Figure [13] lower panel shows the mass ratio of the protostellar outflow 
to the protostar. The ratio increases with the initial cloud mass. As seen in Figure |9l 
there are no significant difference of the mass ejection rate from the circumstellar disk by 
the outflow during main accretion phase among models with different initial masses; in each 
model, ~ 10 — 30 % of the accreting matter is blown away from the circumstellar disk for 
t ^ tff,b- However, the remaining mass of the infalling envelope and outflow driving period 
are different among models. In an initially massive cloud, a large fraction of the cloud mass 
remains in the infalling envelope even after the protostar or circumstellar disk formation. 
The outflow with a wide-opening-angle sweeps and collects the infalling material. Thus, a 
large fraction of the cloud mass swept by the outflow is ejected from the host cloud. The red 
line of Figure [13] also indicates that a massive cloud has a large fraction of the outflowing 
mass. Therefore, the protostellar outflow effectively suppresses the star formation efficiency 
in a massive cloud. As a result, the star formation efficiency in a massive cloud is lower than 
that in a less massive cloud. 



4.3. Star Formation Efficiency 

Figure UM shows the star formation efficiency for each model. In the figure, the diamond 
and square symbols indicate the ratio of the protostellar mass to the initial cloud mass 
(eps = Mps/Mci; O) and the mass of the protostar plus circumstellar disk to the initial 
cloud mass (edisk = [Mps + Mdisk]/^ci, □) at t = to + ts^h, respectively. At this epoch, the 
infalling envelope is already depleted, and gas accretion from the infalling envelope onto the 
circumstellar disk or protostar has already stopped. Thus, in a subsequent evolutionary phase 
of star formation (Class II and Class III phases), the protostar acquires its mass only from the 
circumstellar disk, not the infalling envelope. Thus, solid lines correspond to the lower and 
upper limits of star formation efficiency. The upper limit (square symbol) is realized when 
the entire circumstellar disk finally falls onto the protostar, and the lower limit (diamond 
symbol) is realized when the entire circumstellar disk is blown away or disappears without 
falling onto the protostar by any mechanism such as photo evaporation, jets around the 
protostar, or magnetorotational instability. In fact, the star formation efficiency is expected 
to fall between the lower and upper limits. 

Figure shows that the upper limit of star formation efficiency decreases with the initial 
cloud mass. This is because the mass ejection rate owing to the protostellar outflow increases 
with the host cloud mass. As described in ^ the protostellar outflow is originally driven by 
the first core. Thus, no outflow appears before fi rst core formation in the collapsing cloud. 



The first core has a mass of Mfc ~ 0.01 — 0.1 Mq (jSaigo &: Tomisakall2006l ) at its formation. 
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Thus, in a less massive host cloud, a very slight mass remains as the infalling envelope 
after first core formation. For example, when the initial cloud mass is M^i = 0.05 and 
the first core has a mass of Mfc = 0.04 M©, only 20% [M^nv = 0.01 Mg) of the initial 
cloud mass remains as the infalling envelope. After first core formation, part of the gas 
accreted onto the first core (or the circumstellar disk) is blown away by the outflow. Since 
the outflow is powered by the release of the gravitational energy of the accreting matter, 
no powerful outflow appears unless sufficient accreting matter exists around the driving 
source. In addition, a protostellar outflow with a wide-opening-angle sweeps up the gas of the 
infalling envelope as it propagates into the host cloud. However, when the infalling envelope 
is already depleted, the outflow sweeps up only a slight mass of the infalling envelope and 
ejects it into the interstellar space. The outflow power and the amount of mass swept up by 
the outflow increase with the mass of the infalling envelope. As a result, in a massive cloud, 
a large fraction of the initial cloud mass is ejected from the host cloud and star formation 
efficiency is effectively suppressed by the protostellar outflow. As seen in Table [21 for the 
less massive cloud model N09, 92% of the initial cloud mass accretes onto the circumstellar 
disk and the protostar and only 8% of the initial cloud mass is ejected by the protostellar 
outflow. On the other hand, about half of the initial cloud mass is ejected from the host 
cloud by the protostellar outflow for the massive cloud model N35. 

As seen in Figure [TU the lower limit of star formation efficiency also decreases with 
the initial cloud mass for M^i < 0.26 Mq, whereas it increases slightly for Md > 0.26 Mq. 
The lower limit is determined by the efficiency of mass accretion from the circumstellar 
disk onto the protostar. A less massive disk appears in a less massive cloud and is stable 
against gravity. In such a disk, angular momentum is transferred by magnetic effects such as 
magnetic braking and protostellar outflow and the mass accretes steadily onto the protostar 



(IMachida et al.ll2010bl ). In contrast, a massive circumstellar disk appears in a massive cloud 
and is unstable against gravity. In such a massive circumstellar disk, a non-axisymmetric 
structure appears because of gravitational instability, as shown in Figure [71 This structure 
effectively transfers angular momentum outward, promoting mass accretion from the circum- 
stellar disk onto the protostar. In add ition , the gas accret e s unste adily onto the protostar, as 



described in lVorobyov fc Basu J2006h and iMachida et al.l (l2010al ): these authors pointed out 



the possibility of episodic accretion in such massive disks. Thus, in addition to the magnetic 
effects, the dynamical structure of the circumstellar disk contributes to angular momentum 
transfer in the massive circumstellar disk that forms in an initially massive cloud. As a 
result, the accretion rate from the circumstellar disk onto the protostar increases with the 
initial cloud mass, and thus the protostellar mass and star formation efficiency also increases 
with initial cloud mass. 



Star formation efficiency in the least less massive cloud is e = 0.54 — 0.92 and that in the 
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most massive cloud is e = 0.39 — 0.47. Figure [TH shows that both the lower and upper limit 
of the star formation efficiency tend to decrease as the initial cloud mass increases. This is 
because, in a massive cloud, a protostellar outflow with a wide-opening-angle can sweep up a 
large amount of mass and eject it into the interstellar space. The red line in the lower panel 
of Figure [13] shows that the mass ejection rate owing to the protostellar outflow increases 
with the initial cloud mass. The mass ejection rate owing to the protostellar outflow exceeds 
Mout/Mci ^ 0.3 — 0.5 in a relatively massive cloud. Thus, at most half of the initial cloud 
mass can be ejected by protostellar outflow. 

Since we cannot estimate the mass ratio finally falling onto the protostar from the 
circumstellar disk, the realistic value of star formation efficiency is unclear. However, our 
result indicates that the protostellar outflow contributes greatly to the protostellar mass and 
star formation efficiency in a single cloud core. 



4.4. Initial Cloud Parameters and Spatial Resolution 



In this study, we fixed the initial ratio of magnetic and rotational energies to the gravita- 
tional energy in each cloud. As described in §4.21 star formation efficiency depends strongly 
on outflow properties, which in turn depend on magnetic and rotational energies of the initial 
cloud. Thus, the star formation efficiency described in §4. 21 may be just a lower limit because 
we selected magnetic and rotationa l energies most suitable for driv i ng a powerful outflow 



according to the previous studies (ITomisakal l2002l : iMachida et al. I l2005bl : iMachida et al. 



2008bl ). For example, when a very weak (negligible) magnetic fleld exists in the initial cloud, 
only a weak (or negligible) outflow appears, and higher star formation efficiency is realized 
without significant mass ejection from the host cloud. Even with a slower rotation rate, a 
weak outflow realizes higher star formation efficiency. In addition, the initial strength of the 
magnetic field affects the collimation of the outflo w: outflow appe aring in a weakly magne- 
tized cloud has a relatively narrow opening angle ( ITomisakal 120021 ) . Thus, to investigate the 
relationship between star formation efficiency and protostellar outflow in more detail, we may 
need to investigate cloud evolution in terms of magnetic and rotational energies. However, 
such calculation incurs a very high CPU cost. In this study, we showed that the protostellar 
outflow can suppress star formation efficiency to ~ 25% (~ 50% including the circumstellar 
disk) at most. This indicates that the protostellar outflow greatly affects star formation, and 
we need to consider the protostellar outflow in order to investigate the protostellar mass. 
We will investigate cloud evolution with a large parameter space in the future. 



In this study, to realize long-term evolution of the cloud until the end of the main accre- 
tion phase, we adopted sink treatment instead of resolving the protostar and the structure 
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aro und it (r <^ l^U)- How ever, as described in lTomisakal (120021 ). iBanerjee &: Pudritzl ( 120061 ) . 
and iMachida et al.l ( l2008bl ). another flow component called the high- velocity jet may appea r 
around the protostar. Since the high- velocity jet is well coUimated ( IMachida et al.ll2008bl ). 
it cannot sweep up a large amount of mass in the infalling envelope when it propagates into 
the host cloud. However, about 10% of the accreting matter is expected to be ejected by the 
collimated jet. Thus, the high-velocity jet may further lower star formation efficiency during 
the star formation. However, the mass ejected by the jet may rarely affect the star formation 
efficiency, because a large fraction of the cloud mass is already ejected by the outflow in the 
early stage of the star formation. We need a considerably higher spatial resolution to include 
the effect of high-velocity jets to estimate the star formation efficiency in more detail. 



5. Summary 

In this study, we investigated the impact of protostellar outflow on the star formation 
process. We constructed nine models with different initial cloud masses in the range of 
Mci = 0.015 — 1.5 Mq, and calculated the cloud evolution until the cloud mass is depleted. 
In the calculation, without artiflcially inputting outflow (momentum) to the computational 
domain as seen in any other studies, outflow is naturally or automatically driven by the 
circumstellar disk. As a result of the calculation, we found that a large fraction of the initial 
cloud mass is ejected from the host cloud into the interstellar space by the protostellar 
outflow. Thus, the protostellar outflow signiflcantly affects the star formation process and 
determination of the star formation efficiency. The following results are obtained. 

The protostellar outflow continues to be driven by the circumstellar disk for about 10 
times the freefall timescale of the initial cloud after the cloud begins to collapse. The flnal 
size of the outflow is different among models, because different models have different initial 
central densities and different freefall timescales. In each model, the protostellar outflow 
reaches ~ 500 — 10^ AU far from the protostar in t ~ lOt^^c- When the outflow remains 
inside the host cloud, the outflow extends also to the vertical direction toward the direction 
of travel, because the outflow is anchored by large-scale (host cloud scale) magnetic fleld 
lines. After the outflow penetrates the host cloud and propagates into the interstellar space, 
it extends only in the direction of the travel keeping its width. Thus, in this period, the 
outflow collimation improves to reach ~ 10 — 30. Before disappearing, the outflow has a 
length of 10 times the natal cloud radius and a maximum width comparable to the natal 
cloud radius. Thus, we can expect natal cloud size (or mass) from the length and width of 
observed outflow. 

In this study, since we did not resolve protostar itself, no high-velocity jet appears. 
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However, outflow momentum and momentum flux derived in our simulation well agree with 
observations. This indicates that the outflow driven by the circumstellar disk is responsible 
for total outflow momentum, and high velocity jet rarely affects the star formation efficiency. 
In addition, it is expected that the entrained mass by the high velocity jet can be ignored 
to estimate the mass ejection from the host cloud, because observed outflow momentum can 
be explained only by the disk driven outflow. Although the high velocity jet may promote 
the mass ejection further, a larger fraction of the mass can be ejected only by the outflow. 

The protostellar outflow can eject ~ 10 — 50% of the host cloud mass into the interstellar 
space. The mass ejection rate increases as the initial cloud mass increases; a large fraction of 
the initial cloud mass is ejected in a massive cloud. This is because a massive cloud retains 
a large amount of the infalling matter even after the protostar formation, and outflow can 
sweep and collect a large fraction of the infalling matter when it propagates into the infalling 
envelope. Thus, the protostellar outflow can suppress the star formation efficiency to ^ 50%. 
In addition, a massive circumstellar disk comparable to the protostellar mass remains even 
after the infalling envelope is depleted. Although we need to calculate further evolution of the 
protostellar system to determine the star formation efficiency, the star formation efficiency 
of 30 % may be realized when (a part of) the circumstellar disk is blown away in further 
evolution stage. 

We have benefited greatly from discussions with T. Nakano. This work was supported 
by Grants-in-Aid from MEXT (20540238, 21740136). 
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Table 1: Model parameters and Star Formation Epoch 



Model 


ricfi [cm 3] 


-Rc [AU] 


Mci [Mq] 


Bo [G] 




[s-1 




iff,c [yr] 


to [yr] 


N35 


3 X 10^ 


8700 


1.5 


3.0 X 10- 


-5 


1.3 X 10" 


13 


2.5 X 10^ 


1.0 X 10^ (3.9tff,c) 


N06 


10^ 


4800 


0.8 


5.5 X 10- 


-5 


2.3 X 10" 


13 


1.4 X 10" 


4.5 X 10" (3.2tff,c) 


N36 


3 X 10'^ 


2700 


0.47 


9.6 X 10- 


-5 


4.0 X 10" 


13 


8.0 X 10^ 


3.7 X 10" (4.6tff,c) 


N07 


10^ 


1500 


0.26 


1.7 X 10- 


-4 


7.4 X 10" 


13 


4.4 X lO-"^ 


1.7 X 10'' (3.0iff,c) 


N37 


3 X 10'' 


870 


0.15 


3.0 X 10- 


-4 


1.3 X 10" 


12 


2.5 X lO-"^ 


1.0 X 10" (4.0 iff, c) 


N08 


108 


480 


0.08 


5.5 X 10- 


-4 


2.3 X 10" 


12 


1.4 X 10^ 


4.9 X 10^ (3.5tfr.c) 


N38 


3 X 10^ 


270 


0.047 


9.6 X 10- 


-4 


4.0 X 10" 


12 


8.0 X 10^ 


3.1 X 103 (3.9tff,c) 


N09 


10» 


150 


0.026 


1.7 X 10" 


-3 


7.4 X 10" 


12 


4.4 X 10^ 


2.1 X 10^ (4.9iff,c) 


N39 


3 X 10^ 


87 


0.015 


3.0 X 10- 


-3 


1.3 X 10" 


11 


2.5 X 10^ 





Table 2: Results 



Model 










fps (((liak) 


f mil 






N35 


0.58 


0.13 


0.73 


0.04 


0.39 (0.47) 


0.49 


0.22 


0.03 


N06 


0.23 


0.16 


0.37 


0.06 


0.29 (0.49) 


0.46 


0.70 


0.07 


N36 


0.12 


0.14 


0.18 


0.05 


0.26 (0.54) 


0.38 


1.17 


0.09 


N07 


0.068 


0.09 


0.082 


0.026 


0.26 (0.60) 


0.32 


1.32 


0.10 


N37 


0.043 


0.052 


0.047 


0.007 


0.29 (0.63) 


0.31 


1.21 


0.05 


N08 


0.028 


0.027 


0.024 


0.004 


0.35 (0.68) 


0.30 


0.96 


0.04 


N38 


0.020 


0.019 


0.007 


0.0013 


0.43 (0.82) 


0.15 


0.95 


0.03 


N09 


0.014 


0.01 


0.002 


0.0006 


0.54 (0.92) 


0.08 


0.71 


0.02 


N39 
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Fig. 1. — Time sequence images from the initial state until the end of the main accretion 
phase for model N08. In each panel, the density (color and red contours) and velocity 
(arrows) distribution on the y = plane are plotted with the initial cloud scale. The 
white dashed circle represents initial cloud radius (i.e., the host cloud). The blue line is 
the boundary of the outflow inside which the gas moves outwardly toward the center of the 
cloud (or the protostar) with a supersonic velocity. The elapsed time t in unit of the freefall 
timescale (ts.c) and year is plotted on the upper side of each panel. The white squares in 
each panel denote the outer boundary of the subgrid. 
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Fig. 2. — The magnetic field lines (yellow lines), outflow shape (orange iso- velocity surface) 
and high-density gas region (central red iso-volume) are plotted in three-dimensions. The 
orange surface is the isovelocity surface of Vr = Cs inside which the gas is outflowing from 
the center of the cloud with the supersonic velocity Vr > Cg. The blue sphere corresponds to 
the initial host cloud. The elapsed time t in unit of the freefall timescale (tfT,c) and year is 
plotted on the upper side of each panel. 
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Fig. 3. — Density (color) and velocity (arrows) distributions on the y = plane at t = 
9.39tfr^c for model N06 with different spatial scales. The white dashed circle represents 
initial cloud radius (i.e., BE radius). The red contours denote the outffow inside which the 
gas moves outwardly toward the center of the cloud with a supersonic velocity (thick line) 
and half of the supersonic velocity (thin line). 




Fig. 4. — Density (color) and velocity (arrows) distributions on the y = (left) and z = 
(right) planes at the same epoch as in Fig. [3l The red contour denotes the outflow inside 
which the gas moves outwardly toward the center of the cloud with a supersonic velocity. 
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Fig. 5. — The mass accretion rate (left axis; diamond symbol) and protostellar mass (right 
axis; thick line) are plotted against the time i {—t — to, where to is time at the protostar 
formation) after the protostar formation for each model. The vertical lines in each panel 
corresponds to the freefall timescale of the host cloud at the center (tff^ci solid line) and cloud 
boundary (tff,b) broken line). The initial host cloud mass is plotted by the cross (x) on the 
right axis. 
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Fig. 6. — The mass of protostar, circumstellar disk, outflow and infalling envelope is plotted 
against the time t {= t — to) after the protostar formation for each model. The vertical lines 
in each panel corresponds to the freefall timescale of the host cloud at the center (tfr,c, solid 
line) and cloud boundary (tfr,b, broken line). 
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Fig. 7. — Density (color) and velocity (arrows) distributions on the equatorial plane for 
models N35, NOG, N08 and N09. The elapsed time t in unit of tg^c and year is plotted on the 
upper side of each panel. 
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Fig. 8. — Density (color) and velocity (arrows) distributions on y = plane for models N35, 
NOG, N08 and N09. The elapsed time t in unit of tg c and year is plotted on the upper side 
of each panel. The white dashed circle represents initial cloud radius (i.e., BE radius). The 
white solid line is the boundary of the outflow inside which the gas moves outwardly toward 
the center of the cloud (or the protostar) with a supersonic velocity. The white squares in 
each panel denote the outer boundary of the subgrid. 
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Fig. 9. — The inflowing and outflowing mass ratio of I — /max — 1 and /max — 3 grid against 
the time for model N35, N36, N37 and N38. The vertical lines are the freefall timescale of 
the host cloud at the center {tg^c, solid line), cloud boundary (iff^b) broken line), and 10 times 
the freefall timescale of the host cloud at the center (10tff,c, dotted line). 
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Fig. 10. — The evolution of the outflow length (upper panels) and width (lower panels) in 
unit of AU for all models against the elapsed time normalized by the freefall timescale c- 
The outflow length and width in the right panels are normalized by the initial cloud scale 
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Fig. 11. — The evolution of the outflow coUimation for models N35, N36, N07, N08 and 
N09 against the elapsed time normalized by the freefall timescale of the initial cloud at the 
center. 
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Fig. 12. — The momentum (upper panel) and momentum flux (lower panel) of the proto- 
stellar outflow for all models against the elapsed time after the outflow appears. 




0.1 1.0 

Mcl [Msun] (Host Cloud Mass) 

Fig. 13. — Upper panel: the masses of the protostar (black line), outflow (red line) and cir- 
cumstellar disk (blue line) against the initial cloud mass. Lower panel: the mass ratio of the 
outflow to protostar (Mout/^ps! black line), the circumstellar disk to protostar (Mjisk/Mps, 
blue line), and the outflow to initial cloud (Mout/^ci, red line) against the initial cloud mass. 
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Fig. 14. — The star formation efficiency against the initial cloud mass. The square symbol 
is the mass ratio of the protostar plus circumstellar disk to the initial cloud. The diamond 
symbol is the mass ratio of the protostar to the initial cloud. The star formation efficiency 
in a shadowed area is expected to be realized. 



